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Thirteen targets with mass numbers from 58 to 238 were irradiated with the antiproton beam from 
^q, the Low Energy Antiproton Ring facility at CERN leading to the formation of antiprotonic atoms of 

r—{ ■ these heavy elements. The antiproton capture at the end of atomic cascade results in the production 

of more or less excited residual nuclei. The targets were selected with the criterion that both reaction 
products with one nucleon less than the proton and neutron number of the target are radioactive. 
■ The yield of these radioactive products after stopped-antiproton annihilation was determined using 

gamma-ray spectroscopy techniques. This yield is related to the proton and neutron density in the 
target nucleus at radial distance corresponding to the antiproton annihilation site. The experimental 
data clearly indicate the existence of a neutron rich nuclear periphery, a "neutron halo", strongly 
correlated with the target neutron separation energy B n , and observed for targets with B n < 10 
MeV. For two target nuclei, 106 Cd and 144 Sm, with larger neutron binding energies a proton rich 
' nuclear periphery was observed. Most of the experimental data are in reasonable agreement with 

00 , calculations based on current antiproton-nucleus and pion-nucleus interaction potentials and on 

nuclear densities deduced with the help of the Hartree-Fock-Bogoliubov approach. This approach 
OO ' was, however, unable to account for the 106 Cd and 144 Sm results. 
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PACS numbers: 21.10.Gv, 25.43.+t, 36.10.-k 



I. INTRODUCTION 



There is a growing experimental evidence p|-]l^] that the outer periphery of many stable isotopes of heavy elements 
is composed predominantly of neutrons. Although this is in agreement with simple Jl5[ as well as sophisticated 
p6| nuclear models, its recognition could only rarely be found in the literature [ |l"7||l8| ] . Indeed, the discussions of 
differences between matter and charge distributions were generally limited to the comparison of the corresponding 
mean square radii p9|-[2l|, a quantity much easier accessible to the experiment 22 2j| than the composition of the 



nuclear stratosphere[with the nucleon densities two or three orders of magnitude smaller than the central density). 
Recently, however, the interest in the composition of the outermost nuclear periphery was largely increased when it 
was realised ^6) that the asymptotic behaviour of the nuclear wave function may govern a number of phenomena 
which are expected in experiments with radioactive beams. In particular, a non-homogeneous distribution of the extra- 
neutrons (or extra-protons) in nuclei far from stability may lead to the existence of a marked neutron (proton) halo. 
Such halos would certainly manifest themselves e.g. in low energy transfer reactions induced by radioactive beams 
or, as recently demonstrated p7| , will strongly influence the fusion probability. Before detailed studies of reactions 
induced by neutron rich or neutron poor projectiles are undertaken it is of considerable interest to understand how the 
nuclear periphery of stable nuclei changes its properties as a function of isospin. If this is well described by existing 
theories, one can at least hope that an extrapolation to more exotic nuclei/projectiles (with much larger or smaller 
isospin values than those for stable ones) will be not too far from the reality. 

This example notwithstanding, the composition and extent of the nuclear periphery are evidently of a considerable 
interest by themselves. Governed by the asymptotic behaviour of the nuclear wave function they can be a sensitive 
testing ground for nuclear models going beyond harmonic oscillator boundaries. 

Experimentally, the studies of the nuclear periphery are evidently facilitated if the probes used interact very strongly 
with nuclcons leading to clear signals even from the diluted nuclear stratosphere. Indeed, a substantial part of our 
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present information on the extent and composition of the outer nuclear periphery was obtained via the formation of 
hadronic (pionic, kaonic and antiprotonic) atoms. Two types of experiments can be performed. In the first one the 
X-ray spectra of these "exotic" atoms are investigated p8j-|3l| and in particular the characteristics of the last X-ray 
transition which can be observed before hadron - nucleus interaction prevents further X-ray emission are studied. 
These characteristics are governed by the hadron - nucleus interaction potential depending in turn on the matter 
density where the interaction occurred. In the second type of experiments also involving the formation of hadronic 
atoms, the characteristics of the hadron annihilation products - mesonic or nuclear |L2j - are studied. The 

formation and nature of these products depend not only on the nuclear density distribution but are also related to 
the neutron to proton density ratio. 

We have recently reported fl2] . |l3|| a new experimental study of the nuclear periphery using the formation of an- 
tiprotonic atoms. In this study the nuclear products of antiproton annihilation on a peripheral proton or neutron 
were identified using simple radiochemical methods. It was shown that the outer periphery of the heaviest isotopes of 
all elements studied is composed to a large extent of neutrons. In addition, a strong correlation between the neutron 
to proton density ratio and the neutron binding energy was found. 

In the present paper we give a more detailed description of this experiment together with some new data. 



II. THE METHOD 

After the pure and intense antiproton beam from the LEAR facility at CERN ]32|-|34|] became available the yield 
of the radioactive products formed after stopped antiproton - target interaction could be investigated. Such studies 
were undertaken more than ten years ago |35|-|3^] and were continued thereafter pS-EfJ. Their main objective, similar 
to that for radiochemical work conducted for decades with protons, heavy ions or pions (see e.g. Ref. [[4l]-[l3|) is 
the study of the energy transfer to the target nucleus after the interaction. The unusual character of the antimatter 
projectile, its decay characteristics, the hope for some exotics add to the interest of such investigations. 

Radiochemical experiments with stopped antiprotons were very simple when using the LEAR facility. The extracted 
extremely pure and monoenergetic antiproton beam was traversing a thin scintillator counter (determining the number 
of antiprotons) and was impinging on the target thick enough to stop antiprotons. If the antiproton energy was too 
high a beam energy degrader was used in front of the target. Preferably the target thickness should be larger than 
the antiproton pathlength straggling but this was not always possible with expensive, isotopically separated targets. 

An antiproton entering the target material looses energy by interaction with atomic and conduction electrons. 
When its kinetic energy is comparable to the ionisation energy of the target atoms the antiproton replaces one of 
the atomic electrons and occupies a high-n orbit of this atom. After the capture an electromagnetic cascade develops 
through lower, empty (n,t) states. Spatially, due to the large antiprotonic mass, a good fraction of this cascade occurs 
inside the electronic cloud. The energy released by the cascading antiproton is in the beginning of the cascade (when 
the antiproton is in high n orbits, close to the electron orbits) taken away by Auger electrons. For lower antiprotonic 
orbits the emission of electromagnetic radiation becomes important and this X-ray emission dominates for the lowest 
antiprotonic transitions. The statistical population of the ?-substates at the moment of capture p^-^r| and the 
selection rules for electromagnetic El transitions favour the population of highest-Z levels (=n-l) at the end of the 
cascade. Cascading between these circular orbits the antiproton eventually enters the nuclear periphery where, after 
an encounter with a proton or with a neutron, it annihilates. 

As a result, in more than 95% of the annihilations charged and neutral pions are formed with a multiplicity ranging 
from two to eight and with an average value equal to five. Some of those pions may enter the nuclear volume initiating 
an intra-nuclear cascade [fl7|-[49| and heating the nucleus which subsequently emits neutrons, light charged particles, 
intermediate mass fragments or undergoes fission. The development of these processes evidently critically depends 
on the amount of transferred energy which in turn is related to the number, n; nf , of pions directed towards the inner 
nucleus and interacting inelastically with it. (Due to the strong pion - nucleus potential this interaction is generally 
followed by pion absorption). 

The earliest experimental determination of <nj n t > was presented by Bugg et al. P]. Subsequent experiments are 
discussed and analysed by Cugnon et al. ]5(J and the most recent data are due to Polster et al. As it could be 
expected the average value of iii nt depends on the target mass number A t and <n.i„ t > is about one for A t w 200 and 
less for lighter target nuclei. 

This relatively small value of the average number of pions interacting with the target nucleus immediately indicates 
that the annihilation events with Hi nt = should occur with substantial probability. This is, indeed, confirmed 
by intranuclear-cascade calculations which predict that these "void cascade events" appear in 10-20% of the 
annihilations. After such events cold residual nuclei are formed with a mass equal to the target mass decreased by 
the mass of one nucleon which participated in the annihilation process. Experimentally, such products with mass (A t 
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- 1) were clearly observed with a large yield in the radiochemical studies of the stopped antiproton interaction with 
nuclei US 

In the Nuclear Chart one can find a number of target nuclides for which both neighbouring products with one nucleon 
less than those in the target with mass number A t =N t +Z t are radioactive. For such target nuclei the yield ratio of the 
(N t - 1) to the (Z t - 1) products after the annihilation, easily determined using classical nuclear spectroscopy methods, 
will give information on the neutron to proton "concentration ratio" in the region where annihilation occurred. 

The information on the annihilation site can also be obtained from such a simple radiochemical experiment providing 
the fraction, Y(A t - 1), of annihilation events leading to the (A t - 1) products is determined. This fraction depends 
on the product of the antiproton annihilation probability, W(r), where r is the radial distance from the center of 
the nucleus, with P m i SS (r). The last quantity expresses the "missing probability" jl^] that all pions created during 
the antiproton annihilation miss the inner nucleus, leading to n int = 0. P m i SS (r) evidently increases with the radial 
distance, r. Therefore qualitatively, the larger the Y (A t - 1) value is the farther from the nuclear center the neutron 
to proton ratio is tested. 

The experiment determines the yield ratio of (N t - 1) to (Zt - 1) products and the fraction of events leading to (A t - 
1) nuclei. To interpret the data and to know at which radial distance the method tests the composition of the nuclear 
periphery - the recourse to the theory is necessary. Recent calculations indicate |^3| that almost independently of 
the target mass the antiproton annihilation probability W(r) has its maximum value at distances about 2 fin larger 
than the half-density radius. Its convolution with P m i SS (r) locates the maximum of ni„t = events at a distance even 
about 1 fm larger, making the method sensitive to the composition of the outermost nuclear periphery with nuclear 
densities of about 10~ 2 - 10~ 3 of the central density. We shall come back to these questions at the end of this paper. 

In order to be able to compare the experimental data for any measured target, independently of its neutron to 
proton ratio, we introduce the neutron halo factor defined as 

f perzph = N (P, n ) Im{a p ) Z t 
" N(p,p) Im{a n )N t 

where N(p, n)/N(p, p) is the ratio of produced At - 1 nuclei, and a ra and a p are the p-n and p-p scattering amplitudes, 
respectively [1(1. The ratio Im(a p )/Im(a„) = 1/R np accounts for the ratio of annihilation probabilities. A similar 
halo factor was previously introduced by Bugg et al. ^ who investigated the composition of the nuclear periphery 
determining the characteristics of mesonic products of antiproton annihilation. In his method all annihilation events 
contributed to the halo factor, whereas only the most peripheral events are of importance in the present method, so 
the superscript "periph" is used here. 

Neglecting the corrections mentioned below, the halo factor accounts essentially for the enhancement of the neutron 
or proton concentration over the normal one (i.e. that, which reflects the N/Z target ratio) in the tested nuclear region. 
A halo factor higher than one would indicate an increased neutron concentration ( "neutron halo" ) whereas f/^* p < ~ 
1 would correspond to an enhanced proton density at the nuclear periphery. 

There are, however, corrections which have to be discussed. The above considerations assume that all events in 
which iiint = lead to cold nuclei with mass number (A t - 1). This is not necessarily the case. Microscopically, the 
nuclear periphery may contain some amplitudes of states bound more strongly than the nucleon separation energies 
or fission barrier of these (A t - 1) nuclei. In such a case the antiproton annihilation on these deeply bound states 
would lead to high rearrangement energies pij of the (A t - 1) products, which would subsequently emit particles or 
fission. A simple Thomas-Fermi model estimate [p5| of this effect or more detailed calculations within the Hartree- 
Fock-Bogoliubov model reported in Ref. J53| show that this can affect the halo factor by no more than 30%. 

Another effect which merits consideration is the possibility of the pion charge exchange process which could change 
the final population of (Z t - 1, Nt) and (Z t , N t - I) nuclei. The yield of this process was determined experimentally 
and is discussed in Sec. IVc. 

Ending this section about the experimental method one should mention the article by Bloom et al. [f56| , known to 
us well after the first presentation [ L2|Jl~3| of the method described here. In this article the idea of looking on nuclear 
(rather than mesonic) signals when investigating the composition of the nuclear periphery with exotic atoms was 
already suggested. To the best of our knowledge these suggestions were never realized in the way proposed there. 



III. EXPERIMENTAL TECHNIQUE 



The experimental results which are presented in this paper were gathered during five irradiation series from 1991 
till 1996. The antiproton beam momenta delivered by the LEAR facility are listed in Table |. 

A similar set-up was employed for the experiments with 200 and 310 MeV/c antiprotons. After passing through 
a 18.5 mg/cm 2 beryllium window at the end of the beam tube and after traversing approximately 10 cm of air, the 
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antiprotons encountered the first scintillation counter SI (<f> — 15 mm, thickness 10 mm) with a central hole of 7 mm. 
This counter, used also for the beam focusing, served as an active diaphragm for the antiproton counting during the 
target irradiation. 

A variable thickness plastic beam energy degrader was placed immediately after the SI counter. Its selected 
thickness was calculated to insure that the antiproton beam, after traversing the second scintillation counter S2 (<^> - 7 
mm, thickness 2.5 mm) and a few Al foils (with a total thickness of about 80 mg/cm 2 ), stopped approximately in the 
middle of the irradiated target. Similar Al foils were also placed behind the target. The 24 Na activity produced by 
antiprotons stopped in the backward and forward Al foils was used to determine the fraction of the antiprotons which 
triggered the S2 counter but were not stopped inside the target material. This was due to the fact that the target 
thickness (indicated in Table was generally smaller than the antiproton pathlength straggling (shown in Table | 
and in Fig. |l|) induced during the energy degradation by the moderator, S2 counter and target material. The yield 
for producing 24 Na in stopped antiproton interaction with Al target was determined by independent measurements. 
The average value found was 19.5±1.0 24 Na nuclei produced per 1000 p. The pathlength straggling shown in Table | 
was determined by activation analysis, measuring the activity induced by stopped antiprotons in the target composed 
of a stack of thin foils. 

Depending on the half-life of the (A t - 1) reaction products, each target stack was bombarded by a short (10-15 
min) or long (80-90 min) antiproton "spill", totalizing to about 5xl0 8 - 10 9 particles. Generally more than a half of 
these antiprotons were stopped by the target material, whereas the rest reacted with Al or other target components 
(glue, oxygen). In Table || the targets used in the present work are listed together with the integrated intensities of 
antiprotons stopped in the target nuclei. 

A slight modifications of the above described irradiation conditions were introduced for the last run, where the 
antiprotons with the lowest momentum were available. During this run the SI counter with cj) = 8 mm central hole 
and 1 mm thickness preceded the S2 counter with a 100 /im scintillator. Both these counters were placed in a light- 
tight chamber filled with He gas to decrease both slowing down and scattering of the low-energy beam. Al windows 
of 12 /xm thickness were put at entrance and exit of the He chamber. Eventually, the 6.0 MeV energy delivered by 
the LEAR facility was degraded to 2.8 MeV, with only a rather small energy straggling, as indicated in Table |[ This 
beam was impinging on the target material, with this time large Al foils placed only behind the target to control the 
beam scattering off the targets (diameter 10 mm). 

Shortly after the end of the irradiation (about 2 min in the case of short lived (A 4 - 1) products) the gamma-ray 
counting was started using the HPGe coaxial detectors of about 10 30 % relative efficiency and with an energy 
resolution slightly below 2 keV for the 60 Co 1333 keV transition. In some cases, when the low energy gamma-rays 
were of interest, a planar HPGe detector with an energy resolution of about 600 eV at 122 keV energy was also 
employed. The gamma-ray counting was generally continued for a few days at CERN and later, for many months 
in Warsaw, where in addition a 60% HPGe detector was used. The production yield of (At - 1) nuclei was deduced 
from the absolute intensities of their characteristic gamma lines (followed for a few half-lives) after corrections for the 
decay during and after the irradiation. To this end the antiproton beam intensity as a function of time was always 
carefully monitored. Generally half-lives, energies and branching ratios were taken from the most recent Nuclear Data 
compilations, available via internet connection. In some cases also Ref. ]5^ ] was used. 

For some targets the total number of antiprotons stopped in the target material was additionally determined by 
integrating the mass yield distribution of heavy reaction residues. These distributions were gathered from the yield 
of the radioactive reaction products, as previously described |3^ , |3^| . It was assumed that for non- fissile targets one 
antiproton stopped in the target material produces one heavy reaction residue. The total number of antiprotons 
determined in this way was in general within 10-15% equal to the number obtained using the indications of S2 counter 
with the corrections described above. Figure]^ shows as an example the mass yield distribution obtained for the 130 Te 
target. More information on experimental details and evaluation procedures pertaining to each particular target is 
presented in Refs. p8| , p9[ . 



IV. EXPERIMENTAL RESULTS 



A. Peripheral neutron to proton density ratio 



The final results, used to determine the peripheral neutron to proton density ratios, are presented in Table III. 
These data supersede those previously published by Lubihski et al. |JL|]. The differences are in some cases due to more 
recent information on the decay properties of the (A t -1) nuclei and, in other cases, to the inclusion of results gathered 
during the 1995 and 1996 runs. The second, third and fourth columns of this Table give the absolute production 
yield of nuclei one mass unit lighter than the target mass. The halo factor, , presented in the last column of 
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this Table was obtained from the measured ratio of the produced (N t -1) nuclei to the (Z t -1) ones (column 5) after 
correction for the target Z t /N t ratio and for the ratio of the antiproton annihilation probabilities on a neutron to that 
on a proton, R np . Following Bugg et al. Q the value of the last ratio was taken to be equal to 0.63, in agreement 



with Ref . |6C|1 . Its error is not included in the errors of ff j„ in Table f II 



The data from Table 111 are also shown in Figs. 3-5. Figure || presents the halo factor as a function of the target 
neutron separation energy. The negative correlation, previously observed (T^| for a smaller data sample is confirmed. 
The unusually large error for the 160 Gd target is due to the poorly known absolute transition intensities in the decay 
of (A t -1) products |il|] of this target. The largest value of the halo factor is obtained for the 176 Yb target and is 
discussed below. For two investigated nuclei, 106 Cd and 144 Sm, the measured halo factor is substantially smaller than 
one. As our assumed value of the ratio Im(a n )/Im(a p ) is probably the lowest acceptable (see the discussion in |Q) - 
these results clearly indicate the proton rich atmosphere of these two nuclei. (This would remain true for 144 Sm even 
if one assumes the value R Tlp = 0.48 as obtained in the 4 He experiment ^||). The systematics presented in Fig. || 
shows that the nuclei with smaller neutron binding energy exhibit, on the contrary, a periphery rich in neutrons - a 
"neutron halo" in the terminology introduced more than a quarter of a century ago [^|,^U^]. 

Figure |] shows, as a function of the target mass, another observable of the present experiment, namely the absolute 
yield (per 1000 antiprotons) of the production of nuclei with one nucleon less than the target mass. For all but the 
176 Yb target this yield is close to 10%, without any noticeable dependence on the target mass number. The correlation 
of two observables, the halo factor and the yield of (At-1) nuclei, presented in Fig. || demonstrates again the unusual 
character of the results obtained for 176 Yb. 



B. Isomeric ratios 



If one of the (A t -1) products has an isomeric state with low excitation energy the peripheral antiproton annihilation 
can populate, besides the ground state, also this isomer. Table IV shows the experimentally determined isomeric 
ratios. Not considering the 3~ / 0~ isomeric pair in 152 Eu (where the experimental limit is too high to be significant) 
the high spin / low spin formation ratio is always smaller or substantially smaller than one. 

This is in marked contrast to the isomeric ratios determined for the deep spallation products after stopped antiproton 
annihilation on nuclei |36| , |3"8| ] . In this case, due to the emission of energetic particles in the cascade process, the final 
nuclei acquire a considerable angular momentum, what is reflected by a preferential production of high spin isomers. 
In the case of (A t - 1) nuclei the isomeric ratio depends on the microscopic composition of the nuclear periphery around 
the annihilation site. The contribution of the high spin components is lowered there due to the centrifugal barrier but 
may be increased due to a larger occupation number in comparison with the low spin components. Evidently, in the 
discussion of the isomeric ratio for the (A t - 1) nuclei the feeding of isomers by shorter lived states with excitation 
energies below the particle emission threshold should be also considered. Although such a discussion for all measured 
cases would be outside the scope of this experimental paper, we p resent below in some detail one particularly simple 
example ( 129m,9 Te). We hope that other cases given in Table |y| may be sometimes used to furnish a supplementary 
check for the calculations of the nuclear periphery. 



C. Charge exchange reactions 

The charge exchange process was claimed [f53f as mainly responsible for the larger number of ir~ events in heavy 
than in light nuclei in antiproton annihilation data reported by Bugg et al. (Bugg attributed this difference to the 
neutron halo effect in heavy nuclei.) In the present method the simulation of a neutron halo or at least the increase of 
the observed effect could be due to the transformation of the (N t , Z t - 1) annihilation products to (N t - 1, Z t ) ones via 
7T° — > 7r~ or 7r + — > 7T° charge exchange reactions occurring between annihilation pions and inner nucleus (if they could 
proceed without a substantial excitation of this nucleus). Although such a transformation is evidently undetectable 
by our experimental method, we can estimate its importance looking for the formation of (Zt + 1) nuclei, possible 
only via the same charge exchange processes. 

The determined upper limits for the production of (Z f + 1) nuclei are presented in Table [v|. The comparison of 
the absolute yields of these nuclei with yields presented in the second and third column of Table [ill] clearly indicates 
that the charge exchange effects are generally much smaller than the experimental errors assigned to the production 
yields of nuclei determining the halo factor. Therefore, in agreement with arguments presented in Ref. we will 
in the present work ignore the corrections which could result from these processes. This is not in contradiction with 
37| ], where two (Z t + 1) products were observed. Their absolute yields were, however, not much 
larger than the upper limit values presented in Table M. 



the results of Ref. 1 36 
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V. DISCUSSION 



The method presented in this work has allowed us to correlate in a quantitative way the neutron enhancement 
in the nuclear periphery with the neutron separation energy values, B n . The data presented in Fig. || indicate that 
isotopes having B n smaller than about 10 MeV exhibit a nuclear periphery in which the neutron density is larger than 
that expected from the N/Z ratio of a given nucleus. As the neutron separation energy for the heaviest isotopes of 
all naturally occurring elements is generally below 10 MeV - the enhancement of the nuclear periphery with neutrons 
should be a quite common phenomenon. In the following section we will show that this observation is in a qualitative 
agreement with expectations based on the Hartree-Fock-Bogoliubov approach, a model commonly used to describe 
the properties of the nuclear periphery. To compare the experimental results with theory in a more quantitative way 
the antiproton-nucleus interaction should be considered. This was thoroughly discussed in a previous paper by some 
of us |53| and is also briefly outlined below. 



A. Nuclear periphery from HFB calculations 



In the present paper we concentrate on the description of the nuclear periphery using a self-consistent Hartree-Fock- 
Bogoliubov (HFB) model. We are perfectly aware of the limitations inherent to this model. One of them, namely 
the unability to correctly predict the binding energies, is probably the most troublesome for an approach in which 
the binding energy is used to correlate one of the experimental observables. This was the reason why in our previous 



papers [13 53 we have also investigated another very simple asymptotic density model [15| in which a number of 



phenomenological inputs, including binding energies, was used. However, the agreement with the experiment was 
apparently not improved in comparison with the HFB model. Therefore, at the present stage of this research we 
concentrate on the HFB method. We hope that in future works our experimental data will allow to discriminate 
between different approaches p4| modelling the nuclear periphery. 

We have applied an HFB code which uses the coordinate representation and solves the HFB equation on a spa- 
tial mesh. For the Skyrme interaction, the HFB equation is a differential equation in spatial coordinates. In our 
calculations, we used the SKP version of the Skyrme interaction fl65f . Due to the imposed spherical symmetry, the 
HFB+SKP model is solved separately for each partial wave We used 100 mesh points in the radial coordinate 
in a box of the size of 25 fm. As a boundary condition, we demanded the wave functions to vanish at the far end of 
the box. 

The above described method gives the amplitude of the wave functions which can be followed to large distances from 
the nuclear center. The contribution of each orbital to the nuclear density is given by the square of this wave function 
at a given radius times the occupation number for this orbital. Particle wave functions are normalized in such a way 
that the trace of the hermitian density matrix for protons (neutrons) is equal to the proton (neutron) number. For 
large radial distances the centrifugal barrier suppresses the contribution of high - j orbitals in comparison with those 
with small - j. (This effect is illustrated in Fig. 3-4 of Ref. Jlj}). The centrifugal barrier effect may be compensated 
or overcome by the occupation number for high - j orbitals. As a result, both low- and high-spin orbitals contribute 
to the composition of the nuclear periphery. Figure ^| gives an example of the calculated composition of the nuclear 
periphery for the 130 Te isotope. The antiproton annihilation leads to a hole in one of the indicated orbitals. The 
corresponding hole-excitation state will decay by gamma or conversion electron emission to the 129 Te or 129 Sb ground 
state. In the case of 129 Te the annihilation can also lead to the excitation of the hll/2 isomer. (The 19/2-isomer in 
129 Sb has a three-quasi-particle nature ]66[ ] and, evidently, cannot be excited directly by antiproton annihilation on 
any 130 Te proton orbital. This is in agreement with the experimental result given in Table IV). 

Figure ^| shows only the amplitudes of orbitals contributing with more than 5% to the composition of the nuclear 
periphery. Among states having smaller contributions some are so deeply bound that the antiproton annihilation on 
them leads to excitation energies of (A t -1) nuclei larger than the neutron separation energy (or proton separation 
energy + Coulomb barrier, or fission barrier). These "deep hole" states will evidently decay by particle emission, 
affecting the primary population of the (A t -1) nuclei. The total contribution of these states to the composition of the 
nuclear periphery is also indicated in Fig. ^. Around the antiproton annihilation site (discussed in more detail below) 
this contribution does not exceed 30% for all cases discussed in this work. 

Figure [?] presents the calculated ratio of the (normalized) neutron to proton densities as a function of the radial 
distance for some target isotopes indicated in Table |III| The experimentally observed neutron rich periphery for 
isotopes with neutron separation energy smaller than 10 MeV is qualitatively confirmed by the calculated densities. 
On the contrary, the proton rich nuclear atmosphere observed for 106 Cd and 144 Sm (f/^* p < 1 is not expected by 
the theory). 
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Figure g shows the absolute values of the calculated neutron and proton density for the 96 Zr isotope. Figure |9|, 
finally, displays the neutron to proton density ratio as a function of the radial distance for a series of even Zr isotopes. 
The gradual "build up" of the neutron atmosphere with increasing mass number is clearly shown by this calculation. 
The Figure presents, besides stable isotopes the theoretical expectation for isotopes up to 8 mass units heavier or 
lighter than the stable ones. It illustrates what may be expected when experiments with moderately changed isospin 
will become feasible. 



B. Antiproton-nucleus interaction 



The interactions in question involve three distinctly separate stages. 1) The initial atomic state of the antiproton 
is determined by the long-range Coulomb interaction. At nuclear distances, however, the atomic wave function 
is determined by the centrifugal barrier and an antiproton-nucleus optical potential. 2) The next stage consists 
of antiproton annihilation on a nucleon leading to final mesons. This process lies beyond the present theoretical 
description and must be discussed in a phenomcnological way in terms of the pion multiplicity and energy distribution. 
3) The final pions undergo elastic and inelastic nuclear reactions. The radiochemical method filters these to the elastic 
channel. 

Each stage is fairly complicated and requires an approximate description. Fortunately there are two simplifying 
effects: the strong absorption of the antiproton and pions and the large energy release in the annihilation act. It is 
this energy release that allows to use the closure approximation over the final nuclear states. As a consequence one 
can obtain |53| a semi-classical formula for the partial absorption width on a single nucleon i, 

| $#(Y) | 2 u(Y - X) Ai (X)(l - P dh (X) ) P m iss ( — 7T — )dYdX.. (1) 

1 

Here Y are the antiproton and X are the nucleon coordinates, respectively. The h^n is the reduced mass of NN 
system and pi is the nuclear density for protons or neutrons, given by a sum over single nucleon states, p = 4> a . 
This density is folded over a form factor v(Y — X) that represents the finite extent of the annihilation region (radius 
0.8 fm). The last two terms describe final state effects. P m iss is the probability that final state mesons do not 
interact or interact only elastically. This function excludes the annihilation events which produce highly excited A-2 
or lighter final nuclei. Although calculations of P m iss are fairly involved the net effect is close to a simple geometric 
estimate of an opening angle as seen by pions to avoid nuclear collisions |6^| . Finally Pdh is a model correction for 
the rearrangement energy in the final nucleus. It excludes those captures that lead to A t -1 nuclei excited above the 
particle emission or fission thresholds. 

In the limit P m i SS — 1, Pdh = and pi equal to the total nuclear density Eq. (|l|) produces atomic level widths 
measured in the X-ray experiments. Ima N ^ is believed to represent an effective averaged antiproton - nucleon 
absorptive amplitude and is determined by a fit to atomic level widths and shifts. This number is not needed for 
the studies of neutron halo. What is needed is the ratio R np = Ima n p/ Ima p p which, as indicated above, is taken as 
0.63. (In Ref. J53| the arguments to use the slightly different value R„ p =0.82 were given. In the present paper we 
keep our previous p3[ definition of the halo factor in order to facilitate the comparisons.) For the details of final state 



calculations and the derivation of the intuitively simple formula (1), we refer to Ref. |53|. One problem in the analysis 
of halo factors is to know the atomic state from which nuclear capture takes place. At this stage of research a complete 
answer is not available. One can calculate the distribution of these states for the lower part of the atomic cascade. 
This is possible for states n<20, i.e. states localized between nucleus and electron cloud. In addition, this part of the 
cascade may be checked directly via measurements of X-ray intensities. It was concluded [ p r 3p^ ] that about 95% of 
the nuclear captures take place from levels with only one or two different values of the antiproton angular momentum 
I, although the distribution over n may involve more states. These I values correspond to the circular orbits of the 
so called "upper" and "lower" states of capture. (In the hadronic atom terminology the "lower n" state is the last 
state, which can be observed before the strong interaction prevents further hadron cascading and X-ray emission. 
The "upper" level has the principal quantum number n higher by one in respect to the "lower" level.) The difficulty 
is that these, low n, captures come to at most half of the total nuclear capture rate |39|,f70) . In our calculations the I 
distribution of capture states given by the low n cascade calculations have been used. 

The arguments in favour of this assumption are illustrated in Fig. In this Figure the calculated total absorption 
probability densities in antiprotonic Nd and the corresponding distributions for cold absorption (i.e. those, leading 
to the production of A t -1 nuclei) are shown for several states with different antiproton angular momentum. These 
densities are, for a given value of I, almost n-independent. In the Nd atom the "low ra" capture takes place from 
the 1=6 and 1=7 states. The overlap of the antiproton wave function with the nucleus in these and higher / states is 
localized at the nuclear surface by the centrifugal barrier. On the other hand, for lower values of I another localizing 
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factor arises - the strong absorption of antiprotons. As a result, the radius of maximum absorption reaches some 
limiting value. This penetration blocking indicates that even if some of the antiprotons are absorbed from lower I 
states the spatial scenario of cold capture is not changed significantly. 

C. Comparison of the experimental data with calculations 

1. Annihilation site 

As it was indicated at the beginning of this paper, the experiment determines the neutron to proton density ratio at 
the nuclear periphery, presumably at distances close to the antiproton annihilation site. The calculations mentioned 
above and presented in detail in Ref. [p3| indicate that, almost independently of the target mass, the method is most 
sensitive for densities encountered at distances about 3 fm larger than the half-density radius. The width of the 
annihilation probability distribution (FWHM) for events testing the density ratio is between 2 and 3 fm (compare 
also Fig. 0). 

The experimental observable related to the annihilation site is the production yield of nuclei one mass unit lighter 
than the target mass. This yield, related to the annihilation geometry, depends strongly on the antiproton - nucleus 
interaction, pion - nucleus interaction, pion energy distribution and kinematic correlations. Fortunately, it is rather 
insensitive to the neutron to proton density ratio, which affects mainly the second observable of this experiment. 

The rather good agreement between the experimental data and theoretical estimate for the production yield of 
A t -1 nuclei, presented in Fig. [n], indicates that all factors governing this production are to a large extent understood. 
However, the only large exception, the 176 Yb nucleus, for which the annihilation site seems to be located much farther 
away from the nucleus than predicted by calculations needs further studies (cf. |53]| ). 

2. Neutron to proton density ratio 

Figure ^ shows a comparison between the experimentally determined and calculated values of the halo factor. 
Again, the agreement is fair although slightly less good than for the yield of A t -1 products, indicating the need for 
some improvements in the nuclear density calculations. This is especially true for the two "proton halo" nuclei 106 Cd 
and 144 Sm for which the HFB calculations predict a neutron rather than a proton rich atmosphere at large nuclear 
distances (see Fig. 0). On the contrary, the calculated neutron and proton densities in 176 Yb can be reasonably close 
to reality, the discrepancy between theory and experiment lying mainly in the underestimated annihilation distance. 

3. Microscopic composition of the nuclear periphery from the isomeric ratio 

Lets us discuss, as an example, the case of the 129 Te isomeric ratio. We assume that the decay of the excited 
hole states in 129 Te after the antiproton annihilation on the peripheral neutron proceeds without the parity change 
(El transitions much weaker than Ml and E2 ones). From calculations for the 130 Te target, as indicated above, one 
deduces that the annihilation site is located around r — 8.7 fm in case of the antiproton- neutron interaction. At this 
distance the ratio of the lhll/2 orbital density (main component of the high spin isomer) to the sum of densities of the 
2d3/2 and other positive parity orbitals is 0.26. As no other negative parity orbitals contribute to the nuclear density 
at this distance this is the approximate theoretical expectation of the isomeric ratio, experimentally determined (cf. 
Table 0) to be 0.45±0.15. In a more refined approach one would consider the antiproton annihilation separately for 
each orbital (similarly, as it was discussed in fig]) and calculate the sum of the corresponding negative and positive 
parity contributions. The assumption of the spectroscopic factors S=l should not introduce any appreciable errors, 
as states belonging to the same orbital should be interconnected by fast gamma transitions. 

Similar arguments applied to the 95 Tc isomeric ratio lead to a theoretical value equal to 0.28, to be compared with 
the experimental result of 0.58±0.22. 

VI. SUMMARY AND CONCLUSIONS 

Using the recently proposed new method for the study of the nuclear periphery composition we have investigated 
a number of targets in the mass range from 58 to 238. The experimental data clearly indicate that the periphery of 
nuclei in which the neutron binding energy is smaller than about 10 MeV has more neutrons than would be expected 
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from the target N/Z ratio. As almost all heavy stable isotopes of all elements exhibit neutron binding energy smaller 
than the value quoted above - the neutron rich periphery should be a quite common phenomenon for stable and 
even more for neutron rich radioactive isotopes. However, the observation of this phenomenon may be difficult in 
experiments in which the nuclear periphery is tested with methods less sensitive than offered by strongly interacting 
probes. Indeed, the calculations based on the antiproton - nucleus and pion - nucleus interactions indicate that the 
method used in the present work is sensitive at nuclear distances around 3 fm larger than the nuclear half-density 
radius. At these distances the nuclear density is about 100 to 1000 times smaller than the central one, which leads to 
the evident experimental difficulties. 

In two cases investigated in this work it was found that the outer periphery of the nuclei studied ( 106 Cd and 144 Sm) 
exhibit a proton rich atmosphere. Although this result may be understood intuitively at least for the closed neutron 
shell nucleus 144 Sm, it contradicts the nuclear periphery composition expected from the Hartree-Fock-Bogoliubov 
calculation performed in this work. Therefore it can probably be used as a sensitive test for the nuclear periphery 
models. 

The experimentally determined isomeric ratio for nuclei one mass unit lighter than the target mass will probably 
constitute a similarly sensitive test for modelling the nuclear periphery. If this happens to be true, a number of other 
cases could be investigated in the near future using the new antiproton facility at CERN. 
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FIG. 1. Example of the measured antiproton range distribution in an Al stack, determined by counting the produced 24 Na 
activity. Before impinging on the Al stack the antiproton energy, equal to 21 MeV, was degraded in 800 fim of mylar and in 
2.5 mm of a plastic scintillator. 



FIG. 2. Mass yield distribution for the Te target. The method for the determination of this distribution is given elsewhere 
[40,58]. 



FIG. 3. Neutron halo factor (defined in the text) as a function of the target neutron separation energy, B n . 

FIG. 4. Absolute production yield of isotopes A t -1 having one mass unit less than the target mass as a function of the target 
mass number. 



FIG. 5. Correlation between halo factor and absolute production yield for A t -1 nuclei. 



FIG. 6. Relative contribution of the individual shell model orbitals to the total neutron and proton densities in Te as 
a function of the radial distance, obtained using the Hartree-Fock-Bogoliubov approach. The curve marked Pdh shows the 
contribution of orbitals bound by more than the neutron separation energy. The arrows indicate the most probable distance of 
annihilation on a neutron (upper part) and on a proton (lower part), respectively. 

FIG. 7. Calculated (HFB method) neutron to proton density ratio as a function of the nuclear distance for some investigated 
isotopes. 



FIG. 8. Calculated neutron and proton densities for Zr. 



FIG. 9. The same as Fig. 7 for Zr isotopes. The dashed lines show the density ratio for the unstable Zr isotopes. 

FIG. 10. Total probability densities (full lines) for antiproton absorption from circular orbits in Nd as a function of the 
distance from the nuclear center, together with the corresponding cold-absorption densities (dashed lines). For these curves 
the scale labelled W(i) holds. The numbers at the curves indicate the respective values of the orbital angular momentum. The 
total nuclear density p(r) (left hand scale) is also shown. 



FIG. 11. Comparison of the experimentally determined yield of (A t -1) nuclei with the calculated values. 



FIG. 12. Comparison of the experimentally determined halo factor with the calculated values. 
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TABLE I. Characteristics of the antiproton beam. 



Experiment 


Beam 


Beam 


Measured pathlenght straggling 


Calculated pathlenght 




momentum 


energy 


Target 


FWHM 


straggling 51 




[MeV/c] 


[MeV] 




[mg/cm 2 ] 


FWHM [mg/cm 2 ] 


1 QQ1 


90^ ^ 


99 ^ 




64±4 










197 Au 


56±3 


57.4 


1992 


200.4 


21.4 


27 Al 


22.5±1.0 


21.9 


1993 


200.4 


21.4 


27 Al 


21.9±1.0 


21.9 


1995 


310.1 


49.9 








1996 


106.0 


6.0 


27 Al 


4.6±0.4 


2.9 



a Ref. @. 



TABLE II. Targets and numbers of stopped antiprotons. 



Target 


Enrichment 


Physical form 


Thickness 


p stopped 


Experiment 




[%] 




[mg/cm 2 ] 


in target [10 s ] 




4b Sc 


99.9 


metal 


55.8 


3.4±0.2 


1992 


56 Fe 


99.9 


metal 


43.6 


9.4±0.9 


1993 


58 Ni 


99.9 


metal 


40.1 


7.4±0.3 


1992 


58 Ni 


99.8 


metal 


60.1 


2.9±0.2 


1996 


96 Zr 


85.3 


metal 


35.2 


2.2±0.3 


1992 


96 Zr 


85.3 


metal 


35.2 


5.6±0.3 


1996 


96 Ru 


97.9 


powder + glue 


20.1 


0.7±0.2 


1992 


96 Ru 


97.9 


powder + glue 


20.1 


0.7±0.2 


1996 


nat Cd 


106 Cd (1.26%) 


metal 


37.5 


1.8±0.2 


1993 


nat Cd 


106 Cd (1.26%) 


metal 


39.0 


7.1±0.3 


1996 


106 Cd 


76.5 


metal 


40.0 


2.1±0.3 


1993 


106 Cd 


76.5 


metal 


40.0 


7.6±0.4 


1996 


i2s Te 


98.3 


metal 


353 


9.1±1.0 


1995 


128 Te 


98.3 


metal 


91.1 


3.2±0.2 


1996 




130 Te (33.9%) 


powder + glue 


104 


0.9±0.2 


1992 


130 Te 


99.6 


powder + glue 


37.8 


9.3±2.0 


1993 


130 Te 


99.3 


powder + glue 


93.6 


2.3±0.2 


1996 


144 Sm 


96.5 


metal 


47.2 


3.1±0.2 


1992 


144 Sm 


86.6 


metal 


45.7 


14±3 


1993 


144 Sm 


90.8 


metal 


109 


6.1±0.3 


1996 


148 Nd 


91.6 


powder, oxide 


95 


1.8±0.4 


1993 


148 Nd 


88.6 


powder, oxide 


90 


4.1±0.4 


1996 


nai Eu 




powder + glue, oxide 


105 


3.2±0.4 


1992 


154 Sm 


97.7 


metal 


54.4 


3.2±0.5 


1992 


154 Sm 


98.3 


metal 


107 


1.3±0.1 


1996 


160 Gd 


98.1 


powder + glue, oxide 


40.7 


9.8±2.0 


1993 


160 Gd 


98.1 


powder + glue, oxide 


110 


4.4±0.4 


1996 


nat Yb 


176 Yb (12.7%) 


metal 


102 


0.8±0.1 


1992 


176 Yb 


96.4 


metal 


31.1 


9.8±2.0 


1993 


176 Yb 


96.4 


metal 


114 


8.6±0.5 


1996 


206p b 


95.9 


metal 


85.7 


15±2 


1993 


232 


100 


metal 


100 


23±1 


1991 


232 rp^ 


100 


metal 


4.54 


0.9±0.3 


1993 


232 


100 


metal 


43.0 


1.5±0.3 


1993 


238 tj 


99.8 


metal 


107 


5.2±0.2 


1992 
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TABLE III. Absolute production yield of At-l nuclei, their yield ratio, and peripheral halo factor. 



Produced nuclei 



Target 


Nt—1 

lOOOp 


Z t -l 
lOOOp 


-At-l 
lOOOp 


at i ^i^i ^ 
iv ypn ) 

iv ypp) 


halo 


2^Ni 
Zo ^ 1 1 


45±4 


49±7 


94±6 


0.90±0.12 


1.3±0.2 


40 


111±18 


34±8 


145±11 


3.3±0.6 


3.7±0.6 




39±10 


50±14 


89±16 


0.79±0.17 


1.1±0.2 




33±6 


72±11 


105±8 


0.5±0.1 


0.6±0.1 


JfTe 

bZ 


65±17 


17±3 


82±14 


3.9±1.0 


4.3±1.1 




81±16 


20±4 


101±12 


4.0±0.4 


4.2±0.4 


et 4 Sm 


< 31 


94±20 


110±8 


< 0.4 


< 0.5 


^ 8 Nd 


56±7 


13±4 


69±8 


4.4±0.9 


4.8±0.9 


M 4 Sm 


77±11 


39±8 


116±10 


2.0±0.3 


2.2±0.4 


lt°Gd 


94±23 


17±5 


111±25 


5.5±1.8 


5.8±1.9 




196±28 


26±6 


222±25 


7.6±0.6 


8.0±0.6 


232 rpv 

90 ln 


71±12 


13±2 


84±11 


5.4±0.8 


5.4±0.8 


238 tt 
92 U 


91±7 


19±2 


110±7 


5.8±0.8 


5.8±0.8 



TABLE IV. Isomeric ratios for At-l nuclei. 



Isomer 


Spin 


Encrgy [keV] 


Isomeric ratio 


«Sc(m)/«Sc(g) 


6+/2+ 


271/0 


0.42±0.05 


95 Tc(g)/ 95 Tc(m) 


9 + ; 1 - 

2 12 


0/39 


0.58±0.22 


127 Te(m)/ 127 Te(g) 


11-/3 + 
2 12 


88/0 


< 0.6 


129 Te(m)/ 129 Te(g) 


11-/3 + 
2 12 


105/0 


0.45±0.15 


129 Sb(m)/ 129 Sb(g) 


19 - / 7 + 
2 12 

5(-V0 ( -» 


1851/0 


< 0.02 


150 Eu(m)/ 150 Eu(g) 


42/0 


< 1.3 


152 Eu(2)/ 152 Eu(g) 
152 Eu(g)/ 152 Eu(l) 


(8)"/3- 


148/0 


< 0.016 


3-/0 - 


0/46 


< 5.7 


152 Eu(2)/ 152 Eu(l) 


(8)"/0- 


148/46 


0.11±0.03 


196 Au(m)/ 196 Au(g) 


12-/2" 


595/0 


< 0.02 



TABLE V. Upper limits for the production of Zt + 1 iso- 
topes by pion charge exchange. 



Target 


Isotope 


Spin 


N(Z t +l) 
lOOOp) 


N(A,Z t +l) 
N(A t ,Z t ) 




bb Co 


7/2" 


< 0.24 


< 0.006 a 


96 Zr 


95 Nb(m) 


1/2" 


< 1.3 


< 0.038 




95 Nb(g) 


9/2+ 


< 0.5 


< 0.015 


96 Ru 


95 Rh(m) 


1/2" 


< 6.2 


< 0.12 




95 Rh(g) 


5/2+ 


< 5.5 


< 0.11 


206p b 


205 Bi 


9/2- 


< 0.7 


< 0.009 b 




2 ,)4 Bi 


6+ 


< 0.7 


< 0.009 b 




203 Bi 


9/2" 


< 1.8 


< 0.023 b 



a Production of Z t nuclei assumed to be equal to 40 per 1000 
antiprotons, from systematics based on N(Zt) dependence on 
the neutron binding energy, B n . 

b Production of Z t nuclei assumed to be equal to 80 per 1000 
antiprotons. 
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Range [mg/cm ] 




A 



9 



8 



7 



o 

73 

J3 








7 8 9 10 11 12 13 



B n [MeV] 




A 




r [fm] 




r [fm] 




r [fm] 




r[fm] 




A 




A 



